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ABSTRACT 

This  paper  presents  the  bond  graph  model  of  a  vibrating  ,rotating  shaft  with  disc  mounted  on  it.  Unbalance 
mass  fitted  to  the  disc  cause  the  disc  to  vibrate.  The  bond  Graph  technique  has  been  used  to  model  this  vibrating  system 
and  SYMBOL  Shakti  software  to  simulate  the  model  and  determine  the  frequencies  and  their  corresponding  amplitudes. 
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INTRODU  CTION 

Vibration  in  any  machine  is  undesirable  and  can  lead  to  catastrophic  effects  leading  to  the  loss  of 
manpower,  machine  and  money  and  incurring  a  heavy  loss  to  the  operating  cost  of  the  plant.  If  the  underlying 
causes  are  not  corrected,  the  unwanted  yibration  itself  can  cause  additional  damage.  There  are  several  causes  of 
vibration  which  may  act  alone  or  at  random.  Therefore,  indicating  such  causes  and  bringing  out  their  remedies  is 
very  essential  in  today’s  world.  Keep  in  mind  that  vibration  problems  may  be  caused  by  auxiliary  equipment,  not 
just  the  primary  equipment.  These  are  some  of  the  major  causes  of  vibration.  Vibration  is  a  characteristic  of 
virtually  all  industrial  machines.  On  recording  the  vibration  in  the  machine,  if  it  appears  to  be  above  normal,  then  an 
insignificant  wear  will  appear  initially,  but  it’s  a  prolonged  which  noise  and  heavy  vibration  then  it  calls  for 
immediate  maintenance  action.  Understanding  why  vibration  occurs  and  how  it  manifests  itself  is  a  key  first  step 
toward  preventing  vibration  from  causing  trouble  in  the  production  environment.  Rotating  machines  elements  like 
bearings,  gears,  fans,  rotors,  shafts  are  very  common  in  every  industry.  Failure  of  these  elements  causes  huge 
monitory  losses.  Condition  monitoring  of  these  elements  can  help  in  preventing  the  catastrophic  failure  of  these 
elements,  thereby  saving  downtime  and  monitory  losses.  The  demand  is  increasing  day  by  day  for  increasing  the 
load  carrying  capacity,  enhancing  the  performance  and  service  life  of  mechanical  systems.  Failure  of  machines 
causes  huge  monitory  losses.  Fault  diagnosis  of  mechanical  systems  can  help  in  preventing  the  catastrophic  failure, 
thereby  saving  downtime  and  monitory  losses.  Various  condition  monitoring  techniques  are  available  for  fault 
diagnosis  of  mechanical  systems  e.g.  acoustic  emission,  wear  debris  analysis,  thermo-graph  etc. 
Vibration  monitoring  is  one  of  the  most  successful  techniques  used  for  condition  monitoring  of  the  mechanical 
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systems.  Signal  processing  of  yibration  signals  is  very  important  for  fault  diagnosis  of  machines.  Two  independent 


parameters  simultaneously  direct  the  flow  of  energy  in  and  out  of  any  system  and  surroundings. 

A  bond  graph  is  a  domain  -independent  graphical  description  of  the  dynamic  behavior  of  a  physical  system. 
The  basic  concept  in  the  bond  graph  analysis  is  to  specify  the  flow  of  energy  in  a  system.  In  the  bond  graph  method,  these 
two  parameters  are  defined  by  the  general  terms  of  Effort  (e)  and  Flow  (f). 

The  power  of  the  instantaneous  energy  flow  is  the  product  of  these  two  factors: 

P(t)  =  e  (t).f  (t)  (1) 

Momentum  p  and  Displacement  q  are  defined  as: 

P=  /  e(t).dt  (2) 

q=/  f(t).dt  (3) 

Therefore,  Energy  (E)  from  equation  Eq  (1)  is 

E  =/  e.  f.  dt  =  /  e(t).dq  =  /  f(t)  dp 

Using  the  principal  of  conservation  of  energy,  the  Bond  graph  model  maps  the  flow  of  power  through  a  system. 
By  keeping  track  of  the  power  in  a  system,  the  energy  is  accounted.  Power  is  a  convenient  entity  in  modelling  since  it  can 
be  described  as  the  multiplication  of  two  conjugate  variables  regardless  of  the  engineering  domain  of  its  origin. 
A  bond  graph  maps  the  power  flow  through  a  system  and  simultaneously  describes  the  relationships  between  the  conjugate 
variables  in  each  branch  of  the  system.  In  this  way,  an  accounting  of  all  the  energy  of  a  system,  and  the  conjugate  variable 
relationships  are  used  to  develop  the  describing  equations  of  a  system. 

MATERIALS  AND  METHODOLOGY 

Experimental  Setup  and  Procedure  of  Data  Collection 

An  experimental  set  up  was  built  based  on  model  based  on  the  model-based  technique  for  identification  of 
vibration.  The  experimental  verification  of  the  Unbalance  Identification  of  mass  on  a  shaft  with  a  single  plane  and  two 
eccentricity  locations  has  been  performed  on  a  Machine  Fault  Simulator  Test  Rig  which  was  available  at  CMERI 
Durgapur. 

A  rigid  shaft  considered  to  be  massless  is  mounted  between  two  roller  bearings.  The  distance  between  the  two 
bearings  is  L1+L2,  which  is  60cm.  The  rotor  shaft  and  the  VFD  Drive  motor  are  connected  through  a  flexible  coupling. 
To  measure  the  vibration  in  X-Direction  and  Y-Direction  at  the  two  bearings,  four  accelerometers  are  connected;  two  in 
each  bearing.  The  weight  of  the  disc  is  653  gram  (Ml).  The  disc  position  id  varied  with  respect  to  three  positions  as  shown 
in  Figure  6,  7  and  8  respectively.  Three  unbalance  masses  of  8  gram  (ml),  12  gram  (m2)  and  16  gram  (m3)  are  attached, 
subsequently  to  the  disc  at  eccentricities  6.85cm  (el)  and  6.85cm  (e2)  separately  one  by  one.  Then  the  shaft  is  rotated  at 
rpm  300,  600,  900,  1200  and  1500;  and  simultaneously  the  vibration  readings  (RMS  values)  and  their  phase  values  in  x 
and  y-direction  at  the  two  bearings  were  noted  down.  Artificial  Neural  Networking  Techniques  with  varying  training 
algorithms  have  been  used  to  determine  the  unbalance  plane  of  vibration. 

Mass  Unbalance:  The  unbalance  masses  used  are  8  grams,  12  grams  and  16  grams,  which  are  shown  in  the 
figures  below. 
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Figure  1  Figure  2  Figure  3 


Figure  1:  Unbalance  mass  of  8gram  attached  to  the  rotating  system;  Figure  2:  Unbalance  mass  of  lOgram  attached 
to  the  rotating  system;  Figure  3:  Unbalance  mass  of  8gram  attached  to  the  rotating  system 

Eccentricity 

All  rotors  have  some  eccentricity.  Eccentricity  is  present  when  the  geometrical  center  of  the  rotor  and  the  mass 
center  do  not  coincide  along  their  length.  In  the  present  case,  the  disc  is  considered  absolutely  balanced  where  the 
Geometrical  Center  of  the  disc  coincides  with  the  Centre  of  the  Gravity  of  the  disc.  In  the  current  research  work 
eccentricity  is  the  distance  between  the  Geometrical  Centre  (Centre  of  rotation)  and  the  points  where  the  unbalance  mass  is 
to  be  attached.  There  are  two  locations  at  a  distance  of  6.85cm  from  the  center  of  rotation  (el)  and  at  a  distance  of  4.85cm 
from  the  center  of  rotation  of  the  disc  (e2). 


Figure  4:  The  Rotating  Disc 


Figure  5:  Schematic  Representation  of  Rotating  Disc 


Plane  of  Unbalance 


The  Figure  6,  Figure  7  and  Figure  8  shows  the  different  position  of  the  rotating  disc  placed  on  the  shaft  (left, 
center  and  right  respectively). 


Figure  6  Figure  7  Figure8 


The  design  of  experiments  (Table  2)  is  based  on  the  maximum  possible  combination  of  Unbalance  mass  (m), 
Eccentricity  (e),  Plane  of  unbalance  (p)  and  the  RPM  (r)  values.  The  table  (Table  1)  below  provides  the  maximum  and 
minimum  values  the  independent  parameters/input  to  the  system. 
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Table  1:  Independent/Input  Parameters  and  Their  Levels 


Serial  No 

Independent  /Input  Parameters 

Minimum  Level 

Maximum  Level 

1 

Unbalance  Mass  (m ) 

8gram 

16gram 

2 

Eccentricity  (e) 

4.85cm 

6.85cm 

3 

Plane  of  unbalance  (p) 

15cm  (measured 
from  left  bearing) 

45 cm  (measured  from 
left  bearing) 

4 

Revolutions  per  Minute  (RPM)  (r) 

300 

1500 

Table  2:  Design  of  Experiment 


Set 

Experimental 

Codes 

Unbalance  Mass 

Eccentricity 

Location  of  Unbalance 
Plane  (in  cm) 

(in  Gram) 

(in  cm) 

(From  Left  Bearins  ) 

1 

NFL 

0 

0 

15 

2 

NFC 

0 

0 

30 

3 

NFR 

0 

0 

45 

4 

L8e2 

8 

4.85 

15 

5 

C8e2 

8 

4.85 

30 

6 

R8e2 

8 

4.85 

45 

7 

L8el 

8 

6.85 

15 

8 

C8el 

8 

6.85 

30 

9 

R8el 

8 

6.85 

45 

10 

L12e2 

12 

4.85 

15 

11 

C12e2 

12 

4.85 

30 

12 

R12e2 

12 

4.85 

45 

13 

L12el 

12 

6.85 

15 

14 

C12el 

12 

6.85 

30 

15 

R12el 

12 

6.85 

45 

16 

L16e2 

16 

4.85 

15 

17 

C16e2 

16 

4.85 

30 

18 

R16e2 

16 

4.85 

45 

19 

L16el 

16 

6.85 

15 

20 

C16el 

16 

6.85 

30 

21 

R16el 

16 

6.85 

45 

The  design  of  experiments  as  mentioned  in  the  table  2  has  been  performed  for  300,  600,  900,  1200  and  1500 
RPMvs.  In  total  there  are  105  sets  of  experiments. 

INDICATORS  USED  IN  EXPERIMENT 

NFC:  No  fault  disc  at  Centre,  NFL:  No  fault  disc  at  a  Left  position,  NFR:  No  fault  disc  at  a  Right  position. 

The  number  succeeding  NFC,  NFL,  and  NFR  represents  the  rpm  at  which  the  disc  is  rotating  without  any 
unbalance  mass  attached.  With  unbalance  mass  attached,  the  experimental  code  is  written  in  the  following  format  as 
mentioned  below  one  after  the  other. 

Position  of  Rotating  Disc  (p)  Unbalance  Mass  (m)  Eccentricity  (e)  Revolutions  Per  Minute  (RPM)  (r) 

For  example:C16e  11500  represents  the  rotating  disc  placed  at  center  of  the  shaft  with  mass  unbalance  of  16gram 
placed  at  an  eccentricity  elrotated  at  1500  RPM.  It  is  similarto  the  rest  of  the  experimental  cases. 

Bond  Graph  Model  for  the  Rotating  System  with  Unbalance  Mass 

A  bond  graph  model  has  been  created  for  the  rotating  shaft  -  disc  system  supported  on  two  bearings  which  are 
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The  center  of  mass  of  the  system  is  modelled  by  four  1-junctions  along  the  center  line  of  the  figure. 
Two  of  them  model  the  angular  yelocity  and  two  of  them  model  the  translational  yelocity  with  proper  subscripts  x  and  y 
respectively.  10xjunction  and  10y  junction  represents  the  angular  yelocity  of  the  center  of  mass  of  the  system  about  X  and 
Y  axes  respectively.  Similarly,lx  and  lyrepresent  its  translational  velocity  in  X  and  Y  axes  respectively. 
A  two-port  element  GY  with  modulus  Izz  times  spin  of  the  rotor  models  the  gyroscopic  coupling  between  the  two  one 
junctions  representing  the  angular  velocity  about  X  and  Y  axes.  C  -  R  field  on  the  left  hand  side  and  the  right-hand  side 
represent  the  two  bearings  on  which  the  rotor  is  supported,  where  the  stiffness  and  damping  are  given  by  C  element  and  R 
element  of  the  field.  The  resultant  linear  motion  due  to  angular  and  translational  velocity  associated  with  a  Y-direction  is 
summed  up  in  the  0-junction  appended  to  bond  39  and  34.  TF  element  linked  with  bond  number  44  and  39  converts  the 
angular  motion  about  X — axis  into  linear  motion  in  Y-  direction  corresponding  to  a  mass  center  of  the  bearing  (LHS),  the 
transformer  modulus  being  equal  to  the  linear  distance  between  the  disc  and  the  bearing  center  (Li).  Similarly,  motion  of 
the  center  of  mass  of  the  bearing  associated  with  Y  direction  is  also  modeled.  Same  is  the  case  for  the  right  -hand  side 
(RHS)  bearing  also. 

Unbalance  Mass  at  a  distance  e  (eccentricity)  from  the  center  of  the  disc  is  modeled  as  a  lumped  mass  point  on  the 
disc.  The  motion  in  rotating  y  -direction  is  obtained  by  summing  up  the  linear  motion  from  ly  and  converting  the  angular 
motion  from  10x  by  multiplying  with  Transformer  modulus  ex  at  Oy  junction.  Similarly,  the  motion  in  rotating  x  -direction 
is  obtained  by  summing  up  the  linear  motion  from  lx  and  converting  the  angular  motion  from  10y  by  multiplying  with 
Transformer  modulus  ey  at  0x  junction.  Instantaneous  values  of  ex  and  ey  are  equal  to  e  times  coscotand  e  times  sincot 
respectively.  With  these  0x  and  Oy  junctions  I  elements  represent  the  unbalance  mass.  However,  to  avoid  differential 
causality,  C  and  R  elements  appended  through  1-junctions  on  both  sides  to  model  hypothetical  contact  flexibility  between 
the  lumped  unbalance  mass  and  the  disc. 

Source  of  Effort  (SE  element)  is  connected  to  the  lx  junction  with  a  value  of  -mg  to  designate  the  weight  of  the 
rotor  system.  The  negative  sign  shows  that  the  weight  acts  in  the  negative  x-direction. 
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RESULTS  AND  DISCUSSIONS 

In  this  current  paper,  a  particular  case  of  16  gram  unbalance  placed  on  the  disc  at  eccentricity  el  is  being 
considered.  The  disc  is  placed  at  the  center  of  the  shaft  which  is  rotated  at  1500RPM  using  a  variable  frequency  drive 
motor.  For  unbalance  faults,  frequency  peaks  are  shown  at  IX  RPM.  In  any  rotating  system,  misalignment  and  bearing 
faults  are  evident,  hence  the  Peaks  at  2X  and  higher  frequencies  are  rellected  in  the  Amplitude  spectrum. 
In  the  current  case  frequency  peak  at  25 HZ 

Experimental  Setup  C16E 11500  denotes  that  the  Disc  is  placed  at  Centre  with  16-gram  unbalance  mass  placed  at 
eccentricity  eland  rotated  at  1500  rpm. 

The  cross  Power  spectrum  of  the  particular  signal  as  shown  in  figure  11,  13,  15,  and  17  determines  the  power 
shared  by  a  given  frequency  of  the  two  signals  using  its  squared  module,  and  the  phase  shift  between  the  two  signals  at 
that  frequency  using  its  argument.  It  is  generally  used  for  non-finite  energy  signals  (mostly  not  limited  in  time  signals), 
who  aren't  square-summable.  The  signal's  PSD  is  the  autocorrelation  of  the  signal's  Fourier  Transform,  as  stated  by  the 
Wiener-Khinchin  theorem. 


Left  Bearing  X-Direction  (Channel-2) 


Figure  10:  Channel-2  Amplitude  Spectrum  Figure  11:  Channel-2  -  Cross  Power  Spectrum 


Left  Bearing  Y-  Direction  (Channel-7) 
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Right  Bearing  X-  Direction  (Channel-3) 


Figure  14:  Channel  7-  Amplitude  Spectrum  Figure  15:  Channel  7  -  Cross  Power  Spectrum 


Right  Bearing  Y-Direction  (Channel-8) 


Figure  16:  Channel  8  -  Amplitude  Spectrum  Figure  17:  Channel  8  -  Cross  Power  Spectrum 


Bond  Graph  Results  after  Simulation 

The  bond  graph  model  representing  the  rotating  shaft  on  two  roller  bearings  is  shown  in  Figure  9. 
The  simulated  FFT  results  from  bond  graph  analysis  depicting  frequency  peaks  at  25Hz  is  shown  in  Figure  18-21.  From 
the  figure,  it  is  clearly  eyident  that  the  yibration  in  Y -direction  is  much  more  than  that  of  the  X-direction. 


Figure  18:  FFT  for  Left  Bearing  X-Direction  Figure  19:  FFT  for  Left  Bearing  Y-Direction 
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Figure  20:  FFT  for  Right  Bearing  X-Direction  Figure  21:  FFT  for  Right  Bearing  Y-Direction 

The  values  of  amplitude  spectrum  (from  experiment),  cross  power  spectrum  and  amplitude  values  from  the  bond 
graph  simulation  has  been  tabularized  in  the  Table  3. 


Table  3:  Comparatiye  Analysis  of  Data 


Bearing 

Location 

Frequency 
(IX  RPM)  for 
Unbalance 

Amplitude 
Spectrum  (mm 
RMS) 

Cross  Power 
Spectrum 
(Magnitude 
Linear) 

FFT  Analysis  (Bond 
Graph)  [Amplitude 
(mm)] 

Left  Bearing  X- 

Direction 

(Channel-2) 

25  Hz 

.001201 

.003205 

.002305 

Left  Bearing  Y- 

Direction 

(Channel-3) 

25Hz 

.001101 

.002440 

.001437 

Right  Bearing 

X-Direction 

(Channel-7) 

25Hz 

.001459 

.004201 

.001241 

Right  Bearing 

Y-Direction 

(Channel-8) 

25Hz 

.002478 

.007136 

.001434 

CHANNEL-2  CHANNEL-3  CHANNEL-7  CHANNEL-8 
FREQUENCY  25  HZ  (IX  RPM) 


Figure  22:  Bar  Chart  Showing  Comparatiye  FFT  Analysis 


Amplitude  Spectrum  (mmRMS)  [Esperimental] 


FFT  Analy si  s  (B  on  d  Graph)  [ Amplitud  e  (mm)] 


Cross  Power  Spectrum  (Magnimde  Linear) 
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The  bar  chart  showing  comparative  FFT  analysis  between  the  amplitude  of  yibrations  at  a  IX  RPM,  frequency 
and  the  bond  graph  simulated  output.  The  difference  between  the  experimental  values  and  the  simulated  values  is 
significantly  less  (in  some  cases  it  is  as  less  as.OOlmm)  which  is  clearly  shown  in  Figure  15.  This  validates  the  Bond  Graph 
model  for  the  proposed  rotor  dynamic  system. 

CONCLUSION  AND  FUTURE  SCOPE  OF  WORKS 

In  this  paper,  a  new  approach  has  been  discussed  with  respect  to  the  modeling  of  a  rotating  system  with  unbalance 
fault.  The  Bond  Graph  approach  is  very  convenient  and  can  be  used  for  any  dynamic  system  analysis. 
Using  Fast  Frequency  Transform  in  Bond  Graph  simulation,  frequency  peaks  and  their  amplitude  have  been  determined 
and  compared  with  the  experimental  results.  The  same  technique  can  further  be  extended  for  other  faults  detection  and 
identification  in  a  rotating  system. 
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